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ABSTRACT 

Detailed interpretation of the east-west profile of 
CO-CRUST 1979 refraction/wide angle reflection data has been 
carried out on the basis of horizontal layered model, WHB 
inversion, T’?-X? analysis and synthetic modelling. First 
and secondary arrivals were treated under the constraints of 
Critical points and relative amplitudes to reach a suitable 
crustal model. Two low velocity zones were proposed in the 
uUppernmanaentadlescnust arédepthssofett-14 km and 25-29 km 
which could be related to high conductivity materials in the 
crust. The depth of the Moho is approximately 44.5 km under 
the study area. 

The profile lies entirely within the Churchill province 
and the model derived here appears to be a composite of the 
Swift Current-Suffield model of Chandra and Cumming (1972) 
and the model of Moon and de Landro (1981) which are on 
opposite ends of the profile under cosideration. A 
north-south crustal fault across the eastern section of the 
profile at a longitude somewhat east of 103°W may separate 


the profile area from the regions further east. 
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i INTRODUCTION 

Digececentumyears; “determination “of@crustal velocity 
structure by means of explosion seismology has been a major 
objective of the Consortium for Crustal Reconnaissance Using 
Seismic Techniques (CO-CRUST). During the summer months of 
Io ell emancmmilao ime thrs consortium conducted surveys) in 
southern Manitoba, central and southern Saskatchewan. Near 
vertical reflection and refraction/wide angle reflection 
techniques were employed in these surveys. This study 15 a 
part of this integrated reconnaissance program and uses the 
refraction data for the interpretation of the east-west 
profile in southern Saskatchewan (line-C figure-1.1). Data 
for this profile were collected in 1979 from two different 
shots detonated from Limerick at the western end of the 
profile. The general aim of the study is to establish a 
reliable cross-section of the crustal lithosphere under the 
study area. The following three aspects will be emphasized 
throughout the study: 
(1) establish a velocity depth function and compare it to 
the results obtained in southern Alberta and _ southern 
Manitoba from similar studies. 
(2) look for any seismic expression of the conductive belt 
(see chapter-2) proposed by Camfield and Gough (1975) which 
crosses the same profile. 
(3) detect the so called Riel discontinuity in the lower 
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experiment (after. Kazmierczak, 


In) gohapterna2zy wee discuss background geology and 
geophysics of the study area. A description of data 
acquisition and field procedure has been included in 
ehapter—3.. Chapter-4 describes basic theories used to 
calculate crustal velocity and depth values including their 
upper and lower bounds. Chapter-5 deals with various steps 
involved in the systematic data processing. In chapter-6 
calculated model parameters are interpreted in terms of 
different theories discussed in chapter-4. Finally a 
discussion of the derived crustal model in connection with 
different physical parameters such as TERENOLOG y. 


GOnductiVvity ana porosity etc 1s included in chapter-7. 
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2. GEOLOGY AND GEOPHYSICS OF THE STUDY AREA 


2.1 General Geology: 

The east-west profile of the CO-CRUST 1979 refraction 
experiment is located in the southern part of Saskatchewan. 
The landscape 1S essentially plain and constitutes the 
central part of the greater Interior Platform. Near surface 
sediments of Cretaceous age are essentially flat lying and 
have an average elevation of a few hundred meters above sea 
level. Older and gently rolling landscape has been burried 
Dye OepOSiuS= Clee sqlacial  otil le) sande sand gravel. The 
Phanerozoic sedimentary section is up to 3.0 km thick near 


south-east Saskatchewan. The sediments thin northward to 


the exposed craton of northern Saskatchewan. Thus 
difference is due to varied basin subsidence and 
sedimentation in different regions. Subsidence was 


accentuated at various times throughout the Phanerozoic and 
basin sedimentation was locally affected by movements of 
Structural regions or blocks of the Precambrian basement 
delimited by lineaments as well as by solution of massive 
halite beds of the Middle Devonian Praire Evaporite, which 
also resulted in the collapse of overlying strata. The base 
of the stratigraphic column is a Succession of early 
Paleozoic clastics followed by Ordovician to Mississipian 
carbonate and evaporites and Jurassic to Pleistocene 


craseres sChigure=2 ii)" The clastics of the Cambrian lie 


Hes SoeTapT-ACEt | ‘taued o> aftz Yo sition Ixow- 2609 
snnwartotedewe be # vit qisdsude: s60-98 Berercd et 
ri? esrudtizoes Sk Atele ilistsneas zi squeal 
e7alive aoe) i hate To pians adhaere sit co = fs3an93 


S06 “eibys 3¢ cia (tiadiasaas ope aes i hte tl > 3s 


sae erodn: ets som Rasen ee f, 5 Ani cates. S95 73e%s Te 


ae 
"a 


aslziyvd.nasd 280 ‘sgacz esi Bint fox gts ep tne 196190 ¢ 


site “tevete : bos 582 fe {sissio % erteoge 
seer Waits wa Out “ot qty we aAlindz yIb?isathes ins 
oa ‘Bulswadson wT: zanén (bee s!T .cevaio7eiae2 joss-tsvem! 
: eae sWerias oy id aveddton ie, :-Wegess boeniKe stg © 
éne ‘sonebiedie ae | badis 57 9u5 af oona76! Si Ro 
qav -a>nebtadve -a3i693 ons ct ee hei 
ina store eneet a(t: iuoipqoids, 2smig enh any ts betapas 
20 eansieawo, ris Ber oei34. 2 ds feap! ane Aaizezrsiites aatasal 
serene’ ‘Bi Siig 5 “ii >-70 anabid =o a tone —_ 
Sviewaa le fine 24uelio E pre be teow A aa shame t te. i biel! 
detrte arrogay rants baie ssa nas to chet ota 

obed att + 628138 ohigizeye 9 2 at vi : ooren 
. 268 25 agig@aosue: © bi eaiphaee 

tena mk mie: a 


5 * = a" = ¥ ‘ 
2 : 7 : 
= - 7 7 " : 


Quaternary 
Terfiary 


Mississippian 


Sea level 


-\O000 m 


-2000 m 


Ordovician 


Figure Zed Reference geological column for southern 


-3000 m 


Saskatchewan indicating order of geological systems (arrows 


indicate major unconformities). 


unconformably on the basement. 

Geologic information on the Precambrian comes’ from 
studies of the exposed shield (Bell, 1971) and studies from 
core samples (Burwash and Culbert, 1976). The basement 
under the Phanerozoic appears to be a  south-westward 
continuation of the exposed Precambrian of the shield area. 
Its composition under the Study area 1s basically 
metamorphosed and volcanic-plutonic rocks with weak or no 
CLCMO me Naukeadlyesdurbering Girommmthate Of “other = geologic 
regions in Canada. In fact geophysical data and radiometric 
age determinations indicate that the Canadian Precambrian 
consists of a number of structural provinces (figure-2.2) 
having distinct Gooonveiea! Signatures of their own. 

The Superior province in the exposed shield of eastern 
Manitoba is Kenoran in age and iS dominated by east-west 
trending structures of greenstone and gqneissic belts, both 
reck complexes being extensively intruded by granitic 
Diucons, 

The Churchill province as exposed to the north consists 
mainly of metaSediments and migmatites. Its rocks are of 
Hudsonian age and show variable trend directions. large 
MasseswoL weakly foliavedegranitic rOcKkSware Common. Among 
the regional assemblages near the study area are the Flin 
Flon belt, largely composed of Amisk-Missi 
metasedimentary-metavolcanic units and the Kisseynew belt, a 


Gherssics terrain (figure-2. 3). 
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Geology and Economic Minerals of Canada, 1970). 
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The transition zone from the Churchill to the Superior 
province is known as the Nelson Front and is divided on the 
exposed shield into two subprovinces, the Waboden and the 
Pikwitonei belts. Granulite facies rocks of the Pikwitonei 
and amphibolite facies rocks of the Waboden belt both trend 
northeast-southwest. The Nelson Front is a geologic feature 
o£ tectonic. Significance. It extends at least 1600 km from 
Hudson Bay towards the international border. There is a 
GLeatsdealsOfecontrovensy in sthe lmbernaturesson | both the 
width and location of the boundary zone. Width estimates 
vary from 1 km (Kornik and Maclaren, 1966) to over 200 km 
(Hajnal and Rose, 1979). Most of the proposed hypotheses 
are based on surface geology of the exposed shield and on 
qualitative analysis of potential field data. 

The Wollaston fold belt in north Saskatchewan consists 
of interfingered pelitic schists, paragneisses and granites. 
This zone on the exposed shield has been correlated with a 
zone of very high conductivity known as the NACP anomaly 
(Camfield and Gough, 1977). This long narrow low 
resistivity zone iS traceable across our profile ata 
longitude of about 105°W (figures-1.1 & 2.4). This anomaly 
was related to possible graphite schists associated with a 
fracture zone at the basement. 

Another lineament of some tectonic significance lies to 
PHemnest mbetweene Latitudesm50 N “and =o N \(figure-2.3). 
This is a possible crustal rift discovered by Kanasewich et 


al (1969) which trends approximately east-west. 
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Figure 2.4 Geology of the 


northern part of NACP anomaly 
(after Camfield and Gough, 


nO 7s) 


Further west, the Interior Platform is bordered by the 
Cordilleran orogen (figure-2.2) a geosyncline composed of 
Phanerozoic and Precambrian rocks which has been subject to 
deformations at various times. Folds and faults of the 


Rocky Mountains trend in a north to northwest direction. 


2.2 Gravity and Magnetic information: 

Figure-2.5 (gravity map) and figure-2.6 (magnetic map) 
indicate that our study area lies within a quiet zone of 
magnetic and gravity patterns. The sharp difference between 
western and central regions (at longitude 103°W) may 
represent a rapid change in lithology over a short distance 
probably indicating! a north=South trending fault). Gravity 
and magnetic features of the Superior-Churchill boundary 
zone have been a matter of investigation since 1939. Recent 
studies (Green et al,1979) have shown that there are five 
distinct magnetic zones in the vicinity of our study area 
(figure-2.6). The important features are as follows: 

(1) south central Saskatchewan is characterized by low 
magnetic relief similar to that exhibited by the Kisseynew 
belteot the churchi! 1 sprovincesin;-thegexposed shield. 

(2) southeast Saskatchewan contains a zone of high amplitude 
HMaGrowschomalves Striking north. This zone 1s thought to be 
an extension of the Flin Flon belt beneath the sedimentary 
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(3) along the Saskatchewan-Manitoba border, an area of low 
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Figure 2.5 Bouguer gravity anomaly of southern plains (after 
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Figure 2.6 Magnetic anomaly map covering eastern part of the 


EM 


Study area (after Green et al., 
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magnetic relief occurs, which is interpreted as the southern 
extension of the Churchill-Superior boundary zone. 
(4) in south-western Manitoba, the magnetic anomalies show 
dominant east-west trends, an important characteristic of 
the Superior province. 

THeRSiqnerumeso Methesqravitympatternmshows moresor Less 


Satiilametrend,a1rections and reliets:. 


2.3 Previous seismic studies: 

The central and southern part of the Interior Platform 
has been the location of a number of seismic refraction 
experiments (Steinhart and Meyer, 1961; Weaver, LSG2- 
Maureau, 1964; Hall and Brisbin, 1965; Chandra and Cumming, 
1972: Hall and Hajnal, 1973; Moon and del Landro, 1981) 
Interpreted results show an average crustal thickness of 47 
km in southern Alberta and 39 km in northern Manitoba with 
upper mantle velocities between 8.1 to 8.4 km/sec. The 


Studies also revealed the existence of a sub-basement layer 


On velocity Gr km/sec in the upper crust and an 
intermediate layer of velocity ha km/sec (Riel 
arscontinuity in Canada) in the lower crust. The 


Surtveldaswittmcurcvemt profile (figqure-2.7) which lies. just 
at the western end of our profile will be an important 
reference for comparison. 

liga addition to the refraction studies, detailed 


reflection techniques have been used successfully in the 
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Figure 


southern Saskatchewan and southern 


in Manitoba, 


experiments 


1980). 


(after Kazmierczak, 


Alberta 


area (Cumming and Kanasewich, 1966; Clowes and Kanasewich, 
L2i0pmSeredayei278imerThe profilestieshow al compllexbeearth 
Structure including a number of dipping horizons within the 
eoustwewtthersionificant topographic relief. Combined 
refraction and reflection interpretations revealed the 
existence of a number of major faults in the middle and 
lower enuse which in some cases penetrate the Moho 
discontinuity (Kazmierczak, 1980). The reconnaissance study 
of CO-CRUST 1977 and 1979 led to some general conclusions: 
(ijmthesercust is thicker in the Churchill province than that 
in the Superior province. A major crustal fault might be 
responsible for this difference (Green et al,1979). 

(2) the east-west profiles in southern Alberta, Saskatchewan 
and Manitoba show similar results for the upper crust. 

(3) the Riel discontinuity is evident only in the Churchill 
province. 

(4) the crust thins eastward. 

Recently, combined modelling of seismic and gravity 
data has been carried out by Sprenke (1982). He noted that: 
GigmechesSuperior province 1S characterized by a rather flat 
andueteatureless = Moho “surface sto the north-east, but a 
Strongly tilted surface in =the south-west. 

(2) the Nelson Front is coincident with a fundamental change 
tietheenature of the Moho discontinuity. 
(3) an area of thick crust in the extreme southeast may be 
related to the Williston basin tectonics. 


(4) the NACP zone occupies a gradational zone between two 
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different crustal regions. In the east, crust is of uniform 
thickness while western part consists of a number of 
east-west trending blocks. 

(5) the crust thins towards the east. An anomalous thinning 
of the upper crust occurs in the southern part of the Nelson 
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3. DATA ACQUISITION & INSTRUMENTATION 

The second cooperative seismic refraction survey of the 
CO-=GRUST. group in south-western Manitoba and _ southern 
Saskatchewan was carried out in 1979 and was designed to 
examine the structure of the crust in the Churchill province 
aS a continuation of its previous studies about the nature 
of the Superior-Churchill boundary zone. There was some 
hope of finding out more about a possible major crustal 
fault running in a north-south direction to correlate with a 
prominent feature on the magnetic map (figure-2.6). The 
features which were of additional advantage for conducting a 
Survey in southern Saskatchewan were as follows: 
(1) easy access for shooting and recording due to a complete 
network of available roads. 
(2) unconsolidated surface layer allowing easy drilling. 
(3) a number of bore holes from oil companies’ supplying 
information about the overlying Phanerozoic sedimentary 
section. 
(4) well defined gravity and magnetic expressions of the 
boundary. 

The east-west profile “which is the*recording site of 
Gum davareextends@hrom MimerickinG@lataitudetnao? 38 e4%eceand 
fomgrtuce wa10Gar.¢.4%) MetoseCarlyle (babtitudes49? 327n5' sand 
Wemontude e102: 95-04) Showitacmeline CG ang Ligqure—17i: The 
Umevecsitye ob! Allbertag team Sreconded the main blaststat 
Gundremerbarosss: Melita® and@eRoblin: The latitude and 
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topographic maps. A reading accuracy of 0.5 mm yielded the 
longitude and latitudes with an accuracy of 1.28" (about 25 
meters). Two shots were detonated from Limerick to record 
the east-west profile. Charge sizes were 250 kg and 1375 kg 
respectively. Reciever numbers 18 to 33 recorded signals 
from the smaller charge while seismometers 1 to 17 recorded 
cheslargerGshot.witxsme of e@shots; slatitude)) longitude and 
other shot parameters were compiled by Dr. Z. Hajnal of 
the Department of Geological Sciences of the University of 
Saskatchewan. Shot point to receiver distances were 
calculated using a geodetic computer program by supplying 
the geographical coordinates aS input data. All refraction 
data were recorded on magnetic ies The recorder types 
used were EMR Digital and Analog FM recorders of the 
University of Western Ontario. Each system had its own 
characteristic gain and frequency response over a passband 
of 1-20 Hertz. Each receiver recorded at least 60 seconds 
of data samples after the first arrival. The elevation 
difference between seismometer positions were less than 100 
meters. The recording sites did not deviate more than 4 
kilometers from the linear profile on either side of it. 
The distance between shot point and farthest seismometer 
location was 287.2 kilometers and the receiver spacings 
Hanged meLcOnme/ = tO loa kilometers. All field data were 
Poona terZed gaa bacco Cm OlLem COM sHeLtZe efor the digital 
instruments and between 58 to 61 Hertz for the fluctuating 
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1344 blocks, requiring a tape length of 58.61 feet (6250 
BPI). Each record contained 7200 samples and the first 
three data words provided information about receiver number, 
shot number and orientation of the seismometers 
respectively. The data were written in format (5E16.6) with 
a block length of 1200 samples. Therefore each trace 


consisted of 24 blocks of data on the magnetic tape. 
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4. THEORETICAL CONSIDERATIONS 


4.1 On the use of some basic parameters: 

Applied stress from an explosion results in the 
propagation of the condition of the strain through the 
crustal medium as elastic waves. Depending on the direction 
of movements of the particles of the media there may be 
several kinds of elastic waves. In the case where the 
particles of the medium move in the same direction as the 
direction of the wave propagation, it is called the P-wave. 
There are three main P-wave groups uSually observed in 
seismic crustal studies: Pg, PmP and Pn. These are 
described below. 

COTY POeOrOUD: 

This group represents waves which penetrate into the 
upper part of the basement having a velocity between 5.5 and 
6.5 km/sec. From field examples it has been found that the 
Pg wave may be complex in nature showing different subgroups 
from sedimentary rocks and the crystalline basement. Ina 
continuous profile there may be several segments of travel 
time branches up to the point of intersection with the Pn 
branch and each segment is delayed with respect to _ the 
preceding one. The amplitude of the first phase becomes 
weaker with distance and the next delayed phase _ shows 
Comparatively large amplitudes where it is first observed. 


RepoLeatem distances the later event becomes "a first arrival 
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and the picture is repeated. This splitting feature may not 
Demnorucedm@ fe spacingieoteidetectonseeisaetoo! large. The 
apparent velocity increases with distance and may even show 
a value between 6.5 to 7.0 km/sec. In such cases, these 
waves will have reached the middle or the lower crust. 

(2) PmP group (a reflection): 

These groupt vy bottoms’ “atethelecrustemantie: transitionsor 
boundary with an apparent velocity between 6.8 and 7.6 
km/sec. Pie alemoreveor@lessestrong velocityegradient@orna 
discontinuity exists between the crust and the mantle, a 
reversed travel time branch is generated. This wave appears 
as distinct secondaries and is characterized by larger 
amplitudes near the point of critical reflection. For a 
continuously observed profile a phase correlation for this 
group may be interrupted and split into separated segments. 
Instead of sharp onsets these reflections (or diving waves) 
sometimes may show a long duration sweep (0.4-1.0 sec) which 
indicates a broader transition zone. Hep ieiaolieke, ~ sevstie biates eye 
the PmP may indicate a sandwich like structure of the crust 
(Fuchs, 1969; Clowes and Kanasewich, 1970; Giese, 1976). To 
join these separated segments in a single curve the 
principle of group correlation (see chapter-5) is applied. 
Apparent velocity obtained from the subcritical range is 
usually greater than that Gone Supercritical range. 
Sometimes the velocity at Supercritical range may show an 
apparent value equal to the maximum velocity of the 


Pg-waves. 
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(3) Pn group: 

This wave penetrates the uppermost part of the mantle 
and has ae velocity between 7.6 and 8.5 km/sec. This wave 
group is recorded as a first arrival at larger distances 
(beyond P130-220°Rkms,) « Recently the Pn wave has_ been 
interpreted as a refracted wave penetrating into a medium 
with weak positive velocity gradients rather than as a head 
wave. This branch disappears or is replaced by another 
delayed event after 300-400 kms. The amplitude recorded for 
this wave is often very small and is easily perturbed by 
Local @ivoiser The backward extrapolation of the Pn branch 
usually meets the PmP branch at the point of critical 
reflection. For the cases where Pn branch intersects the 
PmP branch it definitely means that the mantle wave Pn is 
associated with a greater depth than that of the reflected 
PmP and that there exists no plane interface common to both 


Pn and the PmP (Giese, 1976a). 


4,2 Apparent velocity: 

Near surface local inhomogeneities cause small scale 
Pececulani tiecm an Gthemmtraved shimem@curves .@ hormangaverage 
SEruecune mweemonore the variations and ‘compute the least 
EquabesmmeVeOcn LilesimLOre.diffterent@segmentss Whisegivesian 
dveragei: velocity estimation®for the travel time? branch’ under 
consideration. If however, we are interested in the 


velocities corresponding to each point of the travel time 
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curve, an interpolation method (Giese, 1976) can be used. 
Having three data points, the derivative dx/dT is determined 


by a second degree polynomial (Stewart, 1966): 


Cyc Meh a Ro Xe) (Xk) (XK Eee rT) (a) 
PGA GK eX ae) / OX, Hs ) GSS 5) LT EXSY 


Geka ( Xk Gtx, ) WAC SE OS ) CXFSe aT CR) 


Where yeeAe el Se the srdiStance at swhich) the “derivative ~i1s 
bequived va The velocity of interest is’ the reciprocal of the 


slope and T(X) is the travel time at offset distance X. 


4.3 Calculation of Depth: 

The exact determination of the velocity-depth function 
from the travel time data requires a number of 
preconditions. Therefore this task can only be solved by 
introducing some assumptions in order to overcome the _ lack 
Otmeecda ca. Consequently the solution is valid only within 
certain limits. It is useful therefore to have some simple 
methods which allow the determination of limiting values of 
parameters such as the possible maximum and minimum depths: 
(1) Maximum depth: 

From each point of the travel time curve, three 
Quantities can be determined: (a) distance xX (b) 
corresponding travel time T and (c) apparent velocity 1/Va = 


OLAX. Assuming a homogeneous overburden structure of 
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Genstante thickness 82) la point from the reflection branch 


corresponds to a depth 


Ze=ae0.5ktVael/xX — 130° (2) 


For inhomogeneous media, however this value will give the 
maximum possible Z (Giese, 1976). The maximum possible 


overburden velocity can be estimated from 


Mmax savVa t XAT hee (39) 


The maximum depth for progressive branches including head 


waves can be estimated from (Slichter, 1932; Berry, 1971): 


Tmo el Arm) x COS: (Va.T/x) (4) 


(2) Minimum depth: 
For a reversed branch with constant gradient structure 


the minimum depth can be obtained from (Giese, 1976) 


Fe=e0e 5X1 Ol-sind) 7cosd3 sy) 


where, 

6 = sin ‘!{v,/va(p)} 

Vameoevelocity at the surface: 

Va(p) = velocity at the maximum depth of penetration. 


Depths for the points of a progressive branch is calculated 
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from ( Pavlenkova, 1968): 


Ze=e8075X1 (VaeT/x=1) /(VaitT/x+1)} 95 (6) 


where (X,T) are the coordinates of the crossover points 
between the refracted segments. 
(3) Depth from linear segments: 

BOuweneadmewaves i Athe “first ‘arrivals align in an 
approximate straight line segment. A least square estimate 
will give the best apparent velocity and intercept time from 
which the corresponding depths or thicknesses of the layers 
can be obtained. Thickness of the (n-1)th layer of velocity 


Vine leehavingawan. interceprEs time Te(n) 2sequven by .(Dobrin; 


1976) 
FeCl) m=O oT iT) Vin i) Vn vn) Vi ere) 
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This relation can easily be programmed in a computer to 
Solve for thicknesses of multiple horizontal Tayers. 
Alternatively a simple stripping method can be applied. In 


this process a two layer formula (Dobrin, 1976) 


Heer Sr vit) Vi 22) eV OL) aa (8) 


is sufficient. to calculate thicknesses for n multiple 


layers. For the very first head wave branch this equation 


Zi 


essentially gives the-thickness of the surface layer. Next, 
the first intercept on the time axis is used to reduce the 
Bece wveremaicmeche @short!locationeat the base of the surface 
layer. Then using the same formula for the corresponding 
variables of the second head wave branch will give the 
thickness of the second layer. In this way we can proceed 
to determine thicknesses for multiple layers using the 


Simple two layer formula. 


4.4 Exact solution: 

Herglotz (1907) and Bateman (1910) have independently 
given the solution to the problem of inverting travel times 
into velocity depth values. Wiechert (1910) simplified the 
formula and thus the application becomes easier. 
Independently from the final form of the theorem the 
following conditions must be met: 

(1) the function T(X) must be continuous between X=0 and 
S=Koethat 1S, an interruption is forbidden. 
(2) the apparent velocity along the curve, Starting from X=0 
must increase monotonically. 
(3) the derivative daxX/dT must be continuous along the travel 
time curve. - 
There are several different forms to present the solution of 
the integral equation: 

Xp 
ZN pe = (1/m) | cosh (vp/VCH) tax (9) 
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where, 


Vp velocity at maximum depth of penetration. 
Xp = corresponding offset distance. 
V(X) = velocity at arbitrary distance X. 
For low velocity zones an alternate extended form of this 
equation can be used (Gerver and Markushevich, 1967). 
There are several ways to approach the integration of 


the WHB equation for practical purposes. The simplest way 


is to fit the travel time branch to a polynomial 


mere k@ oO eXkt+Oek 4 Cex eet (10) 


The sorder of the polynomial is chosen from the condition of 
best fit. The approximate function .®) is now 
differentiable in any domain of X. Thus the quantities Vp 
and V(X) are easily obtained by simple differentiation of 
Mix eta earstances® Xpe@eand Sex? The integration may be 
performed by simple methods such as_ the Simpson OF 


Trapezoidal rule. 


Ano x method: 

This method is based on the assumption of Straight ray 
Pathseein the overburden and the existence of a first-order 
Guscomtinuity as a refLlecting horizon. For the “subcritical 
DoueronmeolL the retlections branch both conditions are met Lo 
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Straight line with an intercept (22Z/v)? whereas the Slope is 
equal to 1/V*. V specifies the average overburden velocity 
above the discontinuity. 

When applying T?-X? method in crustal studies one has 
to continevoneselfitolrays*refracted as little as ‘possibile 
but which are clearly observed in the records. If all the 
discontinuities are represented in the travel time curves, 
tEheseeybranches Scan be used) toliéstimate™interval velocities 


aeccondingstosthe Dix formublal (Dix, 1955): 


Vi = VEE Te Hae Te eC Toe THe) ae a) 
where, 

Vi = interval velocity. 

Vee Vee=svelocitres® fromeSthe “firstesand®*second™: branches 


respectively. 
tmepeteoe= intercept times®for ™the first and second branches 
respectively. 
Estimations of interval velocities are useful in that they 
can reveal information about velocity decrease or low 
velocity zones within the interval considered, which are not 


readily obtainable from refraction data alone. 
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4.6 Theoretical seismograms: 

Interpretation of seismic data based on travel times or 
apparent velocities alone are not enough. For a complete 
analysis of the data both kinematic and dynamic properties 
of the wavelet must be examined. Synthetic seismograms can 
be used as an aid to the body wave interpretation. Several 
methods have been used recently to compute theoretical 
seismograms for seismic modelling. Examples of the main 
techniques are the reflectivity method (Fuchs and Muller, 
1971), the generalized ray method (Wiggins and Helmberger, 
1974), Disc ray theory (Wiggins and Madrid, 1974; Wiggins, 
1976), Full wave theory (Cormier and Richards, 1977) and the 
WKBJ seismogram (Chapman, 1978; Dey-Sarkar and Chapman, 
1978). Among these several techniques some are 
computationally very expensive and their use is limited toa 
few trial models. We have chosen the Chapman WKBJ 
selismogram because it is cheaper to calculate and is 
effective for limited bandwidth body wave interpretation. 
The solution of a point explosion in a spherically symmetric 
and isotropic medium is obtained by taking the Fourier 
transform with respect to time, of the equations of motion 
and the constitutive equations and expanding the 
displacement and stress parameters in spherical harmonics. 
A ray expansion method is used to solve the radial ordinary 
differential equation with boundary conditions at the free 
Surface, at the center of the sphere and at interfaces 


appropriate to the ray of interest. Then this 
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ray-transformed response is converted to the response at the 
receiver position by evaluating its inverse Fourier 
transform, Only outward travelling waves are considered 
withsa far Sfield “approximation which “requires ‘that the 
property of = athe: media® should change little® within one 
Seismic wave length. 

Practical limitations of seismic interpretations using 
synthetic seismograms are a result of limitations of the 
computational method and the complexities of the observed 
records. For many regions within the earth, the assumptions 
of the simpler theoretical models, particularly of lateral 
homogeneity and no attenuation, are poor approximations. 
Comparison of synthetic with observed records requires a 
knowledge of the source waveform. Moreover, observed 
seismograms require amplitude corrections which may be 
inaccurate because of inadequate knowledge of the instrument 
response and especially the near surface site response. 

Although several dg#fivculties in interpreting 
seismograms remain including insufficient data, ambiguity of 
interpretation and inadequate mathematical treatment, 
sgnificant® insight#to interpretation®eantbe gainedethrough 
analysis of synthetic seismograms and comparison with _ the 


observed field records. 
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4.7 Low velocity zones: 

Gutenberg (1955) suggested the existence of low 
velocieyalayersgin thelerustssiThisshideatwas againypickedwup 
injcherbeginning of the sixties when detailed data.» -from 
explosion studies became available. Landisman et al (1971) 
have given a detailed discussion of the existence of low 
velocity zones in the upper and middle crust. There are two 
possibilities of detecting a velocity inversion within the 
profile under study: 

(1) having a series of wide angle branches in the 
seismogram, average overburden velocities can be calculated 
uSing T?-X? method. Calculation of interval velocities from 
overburden velocities may reveal nearer about a low 
veloctty layensinvthescrust. 

d2)ewaccording# "t0) Snell‘sy law) ta’ /low *velocity«causes an 
TieeSEEUpti Oneofethe traveletime qurve:) AA jump ineethe ftume 
and distance is generated between two branches which is 
spent in passing through the low velocity region. 

By careful study of the observed sSeismogram other 
indirect evidence for low velocity zones may become clear. 
Among these secondary diagnostic features are the strong 
late arrivals after several tens of kilometers (Landisman et 
al ,1971), rapid attenuation of pure head waves with 
distance and absence of associated head waves corresponding 
tOpa particular reflection branch. 

Demonstration of the existence of a low velocity zone 


reduiuesumcLavel Limes from~high quality retraction profiles 
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in combination with reflection times and amplitude data. 
The sialic low velocity zone possibly caused by intrusions 
has been proposed by Landisman et al (1971) to have a lower 
boundary at depths of 10-15 km. Another velocity inversion 
in the middle crust was proposed at 20-28 km depth which 
could be associated with a region of high conductivity in 
the crust (Landisman et al, 1971). Owing to the existence 
Cirelowevelocity=zones the task of inverting T(X) “data into 
the V(Z) domain leads to a nonunique solution. However once 
the time and distance spent in the low velocity layer is 
determined the corresponding portion of the travel time 
branch can be subtracted from the total time and distances 
to re-establish the continuity of the travel time curve as 
demanded by the WHB method. Slichter (1932) and Gerver & 
Markushevich (1966) have given solutions to the low velocity 
problem, Recently synthetic modelling has _ heen found 
Suitable for detecting low velocity layers in the crust 
(Braile and Smith, 1975) when reliable amplitude information 
exists. Bneechapter 6) wes awit (use 7T*+X (analysis .and 
synthetic modelling technique to establish the presence of 


low velocity zones under the study area. 
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5. DIGITAL PROCESSING 


5.1 Presentation of seismic records: 

The simplest method of data presentation is to plot the 
full seismic record using shot instant and location as the 
zero reference. Such presentation requires a large page 
Size, Or a very small scale however, and it is much more 
USeLU ee COmenaVve mon mlargessscale plot) "—retamming "all “the 
necesSary information in it. Therefore a reduced travel 
time plot was used for overall display and processing 


purposes. A reduced time T' is defined as 


gE Geb Aihe (ei) 


where, 

X = offset distance. 

Vr = reduction velocity. 

Vr was chosen to be equal to 6.5 km/sec for most of the data 
DieeiiSe bieShS ae RlnemSignitiucancer ol Vr s1s “that af this 
velocity 1S close to the apparent velocity of the waves of 
interest it will produce a display showing the phases of 6.5 
km/sec velocity arranged in a horizontal line. This kind of 
display facilitates the correlation process (see 
section-5.4). TP oeLllC MVC LOCWiIny MO IMEEITICEDES Ls eCalmenotebe 
specified an average velocity of the crust is used as_ the 


reduction velocity. 
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Theschannel separation should»be as short as possible 
ensuring meleéan®ireadability.. \sUnifermispacingwofethe traces 
was not possible due to variable seismometer distribution. 

No elevation correction was applied for two reasons: 
(1) such correction would require the knowledge of near 
surface velocities in detail. 

Coeehceamounteeotmrcorreation. would be “small since the 
difference between station elevations was small (less than 
100 meters). 

Since we planned to process only the P-wave phases, all 
record sections show only the traces for vertical component 
seismometers. The average spacing between receiver stations 
were about 12 km. A time scale of 0.5 inch/sec was used to 
accomodate the plot within the page size. A time correction 
was applied to display each record from the zero reference 
time, which we selected to be the shot point time. 

In the complete record section, amplitude information 
is a very desirable characteristic to retain but in practice 
it seems to be very difficult due to the following reasons: 
(1) quantitative determination of the amplitude-charge 
Gelatwoneis very difficult to attempt, 

(2) frequency characteristics of individual seismometers are 
different. 

(3) seismometer coupling may be variable. 

Two different methods were tried for normalization: 

(1) correct response of the seismometers were restored by 


divrdimqueach trace by its “gain “factor. Since at large 
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distances the head wave amplitude varies as 1/X? all traces 
had tom, Deemultiplied*by sanvinversetofsthisefactorsrtThenka 
common dividing factor was estimated to bring the maximum 
amplitude of the nearest (to shot) record to unity. 

(2) each trace had maximum amplitude equal to unity. 

In both cases the mean value of each trace was removed from 
the trace before normalization. 

It was evident that even with a variable seismometer 
response, method-1 was more helpful for correlating the 
secondary events. At least it gave some idea about the 
relative amplitudes of the signals. However, a number of 
records required amplifications before approximate relative 
amplitudes could be ensured. This 1S presumably due to 
variations in seismometer coupling, although other causes 


are possible and are discussed later. 


5.2 Spectral analysis: 

AS a consequence of the general availability of large 
computers the technique of spectral analysis of seismic data 
has come into much use during the last few decades. The 
usual seismic record displays the amplitude as a function of 
time (in the time domain) whereas a spectrum shows amplitude 
as a function of frequency (in the frequency domain). The 
phase of a seismic wave can also be displayed as a spectrum 
which is known as the phase spectrum. Evaluation of spectra 


yields very important information about the source and the 
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medium traversed by the signals. 

The application of the Fourier transform converts the 
time domain information to that in the frequency domain. 
For lengthy data usually a moving window is used to estimate 
the Spectra at consecutive window positions. Since 
truncation of data by a window creates side-lobe effects the 
window in the frequency domain is chosen to be as smooth as 
possible. 

Raw spectral estimation was carried out for a series of 
records with inereasingsdistance (figures—-5.1 & 5.2). The 
length of each record taken for spectral estimation was of 
the order of 8 seconds after the first arrivals. Most ot 
the useful seismic information was within this interval. 
Maximum signal energy was found to be restricted to a 
frequency band of 3-15 Hertz. The main peak of the spectrum 
was observed between 3.5-6.0 Hertz. A progressive lowering 
of frequency with distance iS evident. A lower frequency of 
the main band compared to that of the uSually observed 
Signals may be explained as due to the shot hole 
environment. A small change in the frequency content of the 
main peak with distance may reflect a near homogeneous 
structure along the travelling path. There is ae spectral 
beak at 10-12 Hertz jicorresponding to a distance of 58 km 
PROmmenenGniol, POiit mmeunism Deak Gradually dies “out “at 104 
km. Smaller peaks visible in the spectra can be considered 
due to signals generated within the crust by scattering or 


other secondary effects. 
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spectrum for records at distances 12 and 


5.1 Power 


Figure 


31.5 km from the shot. 
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Figure 5.2 Power spectrum for records at distances 59 and 95 


km from the shot. 
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avoeDigutal filtering: 

Filtering means separation, that is, we separate out 
the part of the record we want to study and eliminate other 
Simultaneous but unwanted parts. In order that filtering 
should be successful it is necessary that the wanted and 
unwanted wave motions differ in some property. To improve 
the signal to noise ratio in the seismogram, differences in 
the following properties can be used: frequency spectra, 
apparent velocity, wave polarization and spatial 
correlation. Frequency filtering is a relatively common and 
Simpler technique used in Phuc wCmmOlmLeL Lac tion 
seismology. Although the noise and signal spectra overlap 


over a wide band of frequencies, we selected the frequency 


filtering because of its SimplLic yey and reasonable 
effectiveness. zero phase shift eight pole recursive 
Butterworth & Bessel filters were used BOE bandpass 


filtering (figure-5.3). A discussion about the comparative 
effectiveness of these filters can be found in Chiu (1982). 
It was shown that: 

GP=the=Besselfvitem has very Vittle “distortion and time 
shift effect on the first breaks but does not suppress noise 
as efficiently as the Butterworth filter with the same 
general characteristics. 

(2) the Butterworth filter can suppress the unwanted parts 
more efficiently than other similar filters due to the rapid 
falloff in response outSide the acceptance band but it 


Gistores the first arrivals significantly. 
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Figure 5.3 Example showing frequency filtering using an 


eight pole Butterworth filter. 
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We therefore used the Bessel and Butterworth filters 
respectively to study the first arrivals and the secondary 
events in the seismogram. THe wepurpose Ofathesril tering 
process was two fold: 

(1) to get rid of the ground noise without affecting the 
first breaks and thus get a seismogram for overall display 
purpose (figure-5.4) 

(2) to enhance the secondary arrivals (figure-5.5). 

For the first case a Bessel filter with a band width between 
4 to 15 Hertz was found suitable. For the second case a 
6-13 Hertz Butterworth filter was used. The later process 
improved the secondary arrivals significantly and aided in 


correlating the later events. 


5.4 Correlation: 

Correlation implies fitting curves or straight lines to 
the same event on different traces. When correlating, the 
seismologist looks for three important properties: - (1) 
character (2) amplitude and (3) frequency of the waves. 
Phase correlation 1s difficult 1£ recording spreads are not 
continuous. Geophone spacings of less than 1 km is 
reasonable for this purpose. But if however this condition 
is not met, one has to rely on group correlation techniques, 
although velocity estimation from this method may differ 
from that of the phase correlation (Giese, 1976). 
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Figure 5.4 Vertical component seismogram of line-C from 


CO-CRUST 1979 refraction experiment (Filter 4-15 Hz Bessel). 
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Figure 5.5 Seismogram of line-C showing correlated 
reflection events from CO-CRUST 1979 refraction data (Filter 


6-13 Hertz Butterworth). 
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Pediat vOiTmeveloecityesis Yclose = toethe velocity of interest. 
The choice of channel separation is a compromise between 
clear readability and smallest possible detector spacing. 
Usually the first peaks or troughs after the beginning of a 
Particular event are correlated. Then a small correction of 
the order of 0.10 sec is applied before calculating depths. 
This correction is also applicable to the intercept times. 
Thema coulertametOr ethe “correlation “process sare vas 
PoOliowsis 

(1) amplitude of events being correlated must exceed those 
of the noise. 

(2) travel time branches must be of some length if they are 
to be identified with certainty. | 

(3) the apparent velocity must show values within a possible 
and reasonable range. 

One has to be very careful when correlating events, 
because once the correlation has been made many aspects of 
the final model are fixed. This is eSpecially true for a 
Single profile interpreted in terms of a plane horizontally 
layered model, since correlation here exactly defines the 
velocities and ultimately the layer thicknesses, except in 
Cases Ot ambiguity of vvlow “velocity zones; Amplitude 
information is also important. If relative amplitudes are 
not reflected exactly, one might correlate wrong events 
together. Attenuation causes progressive lowering of 
frequencies. Therefore when correlating, it should be 


remembered that the wave shape might change with distance as 
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well. The fact that we have used two different shots for 
recording a single profile might cause amplitude differences 


due to the variation in the charge sizes. 
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6. RESULTS AND INTERPRETATION 

[he interpretation® of 9 thé refractuon S@datasfrom™the 
east-west profile includes the analysis of results obtained 
from the horizontal layer model, WHB inversion, T’?-X* method 
applied “to possible reflected phases and finally the 
Synthetic modelling. The final model is a combination of 
all the separate models which best fits the observations. 
Although it 1s not possible to explain all the features of 
the observed seismograms, an attempt was made to accomodate 
the maximum possible numbers of important features in the 


final model. 


6.1 Horizontal layer model: 

The direct arrival can be identified only on the first 
record at, 11-./55 km from the shot point. Due to the large 
spacing between shot point and the first seismometer 
location, we probably missed detail and convincing 
information about the direct wave. From the slope of the 
hineejOtning the only direct: phase to the origin) of ‘the 
plot, we estimated the sedimentary velocity to be 3.64 
km/sec. We compared this value to the available well data 
obtained from nearby locations. The formations above the 
basement indicated an average velocity of 3.2 km/sec for the 
sedimentary layer. The difference is small but we preferred 
the value 3.2 km/sec for further analysis because velocity 


estimation from a two point segment is rather unreliable. 
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The next arrival is more well defined. Events were 
recorded clearly up to 70 km and permit an exact phase 
Gaynelaticnzaelhesvelocit yepicked@upetrom the arrivals upato 
Soeckmnevy seid Ean average value of 5.97 km/sec. The rest of 
the arrivals indicate a slightly increased velocity, the 
least Square estimate of which is 6.2 km/sec. Titeeics 
possible that these two sets of arrivals correspond to a 
Single basement layer with a positive velocity gradient in 
the lower part, but definite conclusion about this was not 
possible due to sparse receiver locations. The first 
arrivals after 70 km are very weak and delayed by about 0.1 
second. Careful study of these arrivals were very important 
because of the following reasons: 

(1) ambiguity of interpreting these arrivals as a 
continuation of the 6.2 km/sec phase due to the observed 
delay. 

V2 munerweaksphases could also be interpreted aS a part of 
the wveflectrortubranch® ihrompiatpossiblesinteriace ate i tekm 
depth. 

(3) since the NACP anomaly crosses our profile at about 
105°W longitude, a study of the signal amplitudes between 50 
ad ee CO terms ebrom GhesSshobmpoantewaseparticularly amportant 
Homvemttymtheapossibleweitectstcoti athe iconductien’ zoneweon 
these arrivals. 

We therefore took special precautions to read the exact 
Sieimeucdes multi ime 00 Ski lomelecce stromemthees Shoue@erhree 


different normalizing methods were tried: 
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(1) each trace was divided by its gain factor to restore the 
exact response. The effect of 1/X*? attenuation of the head 
waves at large distances was minimized by multiplying each 
trace by X*. Each trace was then normalized with respect to 
the average ground noise amplitude observed before the first 
arrivals (figure-6.1). 

(2) this method was similar to the first one except that the 
traces were normalized with respect to the average noise 
level obtained from the far end section of the traces, on 
the assumption that shot generated noise should be of equal 
amplitude on each trace. 

(3) this method tried to minimize the variable seismometer 
responses between 70-100 km from the shot point. Amplitudes 
of the refracted/reflected arrivals from the next layer were 
adjusted to a constant value to observe the resultant 
variations on BiGw Bi eUSt aria valesSldqmals 6 Of “<intenest 
ChEoure- 6.2). 

Detailed study of the amplitudes revealed that the weak 
pulses are reversed in phase with respect to the 6.2 km/sec 
arrivals and can be interpreted as reflection arrivals from 
Armntertace ae fi kine (Obtarned from T?=x* method discussed 
Wmesect Ton 6.3 eC COL di nGmtomsratlemend Smitha Sm5) this 
however indicates a low velocity layer in the lower basement 
(figure-6.2 & Grsa)n Therefore the horizontal layer 
PHECEDLreLaLion was no lOongem Mvalid eior “this “region, We 
therefore used the T’-X? analysis and synthetic modelling 
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Figure 6.1 Normalization of traces with respect 


average noise level before the first arrivals. 


(KM) 


De Sel te 


to 


50 


the 


te 


f1xs8 23i3,)0d coie ted 


fis 
(KM) 


30 


DISTANCE 


7 


ao 


=) 


O 
Oe is 


| OSIS oe Srv Sea iisy Ident 


Figure 6.2 Normalization of traces keeping amplitudes of the 


prominent phases constant. 
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Figure 6.3 Synthetic showing arrivals (PLR) from top of the 


low velocity layer in the basement. PC, PCR & PCR1 specify 


the refracted, reflected and multiple reflected branches 


from the sub-basement layer (after Braile & Smith, 
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6.4. The depth of the basement computed from the average 
velocity of 6.1 km/sec layer and the observed intercept time 
(using equation-7, chapter-4) corresponds to 3.03 km. The 
uniqueness of this value can be justified from the well log 
information as well. We computed theoretical intercept 
times from the total time spent in each layer and compared 
it to the observed intercept time in the seismogram. The 
two values were 1.45 and 1.34 seconds as determined from 
well log data, and as observed on the seismic profile. Tet 


Werealiow wav correlation eprom cofs0ell second (which also 


includes the effect of a possible gradient zone), this 
implies a fairly good agreement. We list the computed 
parameters from well logs and the observed data in 


table-6.1. Most of the well logs were situated northward, a 
few tens of kilometers away from the profile. Since the 
Williston basin gets shallower to the north, a shallower 
basement to the north of our profile was confirmed by the 
well data. The increased depth of the basement as observed 
southward by seismic studies therefore does not contradict 
the well data obtained from the north. 

A secondary phase corresponding to an apparent velocity 
of 4.1 km/sec was identified in the observed seismogram. 
(figure-6.4). Since this phase is absent in the radial 
component seismogram, it could be interpreted as a P-wave 
arrival from the lower sedimentary section. The computed 
thickness and depth of this layer are 2:0 and 1.5 kms 


respectively. Since we lacked further evidence to confirm 
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Figure 6.4 Correlated secondary arrivals corresponding to a 
possible layer of velocity 4.1 km/sec (filter 3.5-8.5 Hertz 


Butterworth). 
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this layer we ignored this secondary event while calculating 
the final model. 

The next set of aligned phases are very prominent from 
70 km to 140 km. The apparent refractor velocity from least 
Square analysis is 6.49 km/sec. These arrivals are followed 
by another onset having an apparent velocity of 6.7 km/sec 
which appear as first arrivals between 150 and 230 km from 
the shot point with a comparatively weaker energy content. 
The refraction arrivals from the Moho region start to appear 
at 225 km from the shot point. The wave pattern has a very 
low frequency content and is perturbed severely by the 
ground noise. However, we were able to pick up the Pn 
arrivals from four geophone locations. The corresponding 
least square velocity is 8.12 km/sec. Continuation of this 
phase towards the corresponding reflection branch (PmP) 
could not be traced possibly because of its weak character. 
Considering the first arrivals from all the evident 
refractors we estimated their depths from the respective 
intercept times uSing either equation-7 or equation-8 thus 
deme rngrcurerirstes ough model (model) Mi, stigure—-6.5) from 
the horizontal layer interpretation. We also computed the 
maximum and minimum bounds of depths using equations-4 & 6 
in Gnapeer—4. The computed results are included in 
table-6.2. From the table we see that the total thickness 
of the crust is 44.5 km under the study area. The estimated 
horizontal layer model lies within the limits of the maximum 


and minimum values. Apart from the observed discrepancy 
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TABLE:6.1 Comparison of results from well logs 
and seismic observation. 


parameters well data selsmogram 
(1) 2-way intercept 1.45 sec 1.34 sec 
time 
(2) average velocity 3 Zia kin, sec 3.64 km/sec 


of sediments 
(3) depth of basement 2.66 km 370 skim 
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TABLE:6.2 Nominal parameters as calculated 
for the Horizontal layer model. 


Layer Velocity Thickness Depth Min-Max 

km/sec km km depth km 
1 Ch 2a 0.0 On Oa ONG 
2 6.0 Piro) Fa eys) Poy aime iat 
3 Gree Sha 2 oll) Sie gs ae 
4 Goo 5.4 14.9 ee om 
S Gree? 24°, 2 AW, &) he a2 aS 
6 sea oa 44,5 4076-540 
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Figure 6.5 Model M1 from horizontal layer interpretation. 
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between theoretical and observed critical points of the PmP 
branch (132 and 152 kms respectively), the results from the 
horizontal layer model can be considered fairly convincing 


UpelOmeNnLSepolnt. 


6.2 WHB inversion: 

The plane layered interpretation from the first arrival 
data assumed constant velocities within each layer. This 
assumption is only an approximation to the more general case 
where velocity usually increases with depth due to 
overburden pressure (Birch, 1964). Therefore if we assume 
that the horizontal layer model (considering straight ray 
paths) gives only a minimum estimate of the depths, we must 
introduce another method which will give exact maximum 
depths. This method is the WHB inversion and is equivalent 
to the method of getting approximate upper limits of depths 
discussed in chapter-4. 

All the first break times across the seismic section 
were approximated by a third order polynomial curve fitted 
in the least squares sense. All data were reduced to the 
base of the sedimentary layer before fitting. Then we 
differentiated this approximate function to estimate dT/dx 
GomLeC ll peCCaDymUNCEELeGtoumReCUdELONSJe ine .chapter-—4™ sand 
computed the corresponding depths. The integration was 
performed by using the Trapezoidal rule. The results 


obtained from this method were similar to the horizontal 
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Figure 6.6 Model M2 from WHB inversion. 
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WayeremodelWexceprwior the regiontof, thevlower",cruset!. The 
depth of Moho obtained from the WHB method was 52 km whereas 
the one computed from horizontal layer interpretation was 
£2 oumkme (tables 6.2) % The difference is presumably due to 
the effect of the gradual increase of velocity with depth in 
the WHB method. Therefore it is clear that the introduction 
of-a velocity gradient in the lower crust will increase the 
depth of the Moho to an amount not exceeding 52 kilometers. 
The model M2 as derived from the WHB inversion is shown in 


PugubesomG. 


6.3 Analysis of reflection branches: 

One of the principal aims of this experiment was to 
record wide angle reflections. Indeed we recorded very 
prominent reflection arrivals. Both Ssubera tical and 
Supercritical branches were observed in the seismogram. The 
Simplest T*-X* method was used to analyse these arrivals. 
At the beginning of the processing work, the reflection 
branches were partly obscured due to high noise level even 
aipenmea 4-15 §Hertz bandpass "filtering. We decided to 
Suppress ene lows frequencies sand "give preference to "the 
comparatively high frequency reflections. With a band width 
of 6-13 Hertz the reflection arrivals became fairly easy to 
ebserve (figure-o.0, chapter—5) 9A group correlation method 
VoSeoDD UNCC eUOnDIC Wm UNCmner LeCUlOnm™Dranches. Rats MUR ee 
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Gorresponding —depths® sat 9 11,5 1327,5 19.5; 24, 27 ande43.4 
kilometers (table-6.3). The interval velocities were 
computed from the Dix formula (equation-11, Chapter-4). 
BKeltowing Greses(1976) a correction of 910% & 3% for the 
computed depths and velocities were applied for the wide 
angle branches. Since correlation of secondary events 
involve the risk of large errors we did not intend to 
include these computed values in our final model. Rather we 
used PheMenadSeeaneeald. tO "a ssdUalitative sinterpretation 
justifying the presence of the reflectors in the crust. The 
EGbLECCOES seat. )13./, 019.5 anduts.4 kilometersacorrespond to 
the refractors mentioned in section-6.2 at similar depths as 
determined by the horizontal layer model (see table-6.2). 
Tiesteriecvom a. 24° km “ofmmdepth ~does: Bnomlsappean aS “a 
Eebractor. This’ i1s= due to ga velocity sinversion at this 
depth (see computed interval velocities in table-6.3) and 
hence would not be observed in the head wave data. The 
interface at 27 km however should refract Seismic rays, 
Since it apparently represents an increase in velocity as 
determined by interval velocity calculation. There iS no 
elear explanation for the apparent absence of this refracted 
phase in the observed seismogram. Moreover the observed 
reflection branch corresponding to 24 km depth was so 
Promimlent thac we COULGMnOLe Ignore it -or takes it “as) a 
multiple reflection. ihe is possible that the _ weak 
refraction arrivals from the 27 km interface were obscured 


due to the low signal to noise ratio beyond 200 km. Another 
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TABLE:6.3 Results from T?-X? analysis 


Layer Interval velocity Depth Min-Max 

km/sec km depth km 
1 DAS 00.0 OO 0=0 010 
2 6.4 SS rem 2 ee Oaal Ore, 
3 Ge6 We S 19. 522470 
4 Spe 24.0 RBs I= Rahm f) 
5 Grad 2h. 4 Le Os 2 eo 
6 oie | 43.4 42.0-49.0 
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possible explanation requires a gradient zone between depths 
of 24 and 43.4 kilometers so that a head wave from 27 kms of 
depth is no longer possible, instead a retrogade branch is 
generated by the corresponding diving waves. For the later 
case the layered model in the lower crust will be only an 
approximation to the real situation. The crustal model 
obtained from the T?-X? method (model M3) is plotted in 
figure-6.7. The average velocity of the crust (6.39 km/sec) 
obtained from this method is in good agreement with previous 
results (Sereda, 1978; p.193). We also computed the upper 
and lower bounds of the depth values for the reflection 
branches uSing equations-2 & 5 from chapter-4. The result 
PSs neludced Sinm table-G.3. silts interesting» to mote that 
the calculated minimum values are almost equal to the T?-X? 
results. This implies that these two methods are equivalent 
BOmeCAChMOLNeC Hmm USnalSOmMeanSmenacmtne lx —Methnod yw ild 
always give a minimum eStimation of the depth values. 
Although the computed depth and interval velocities were 
used as a diagnostic means to Support the low velocity and 
the reflection branches, the exact values will be fixed by 


the synthetic modelling technique described below. 


6.4 Synthetic modelling: 
Interpretation waS carried out by model calculation 
using the synthetic seismogram formulation described by 
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Figure 6.7 Model M3 from T?~X? analysis. 
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inversion can be used as an input for the synthetic 
seismogram. Usually the synthetic obtained from the first 
model does not fit every feature of the observed data. 
Therefore one then tries to fit between the calculated and 
the observed. travel times by changing the first model in a 
SystematbicGway ai MBy wea Gtrialeiand Gerroretprocessuiones can 
eventually reach a good agreement. We did not have a sure 
Strategy for changing the model ina certain direction to 
reach better and better results, but we did find that by 
applying few changes at once and starting from the surface 
downwards one obtains results most quickly. In any case 
uniqueness will alwayS remain questionable. By using a 
trial and error method in changing the models one finds that 
enesustalwaysi pushed@asin¥ia Tcertam: tdirectiony Se that las, 
independent of the type of model one starts with, one will 
end up with a particular class of model. This however does 
not imply that the final model will truely represent the 
aetua bWcenditions mrt WwElitistailebe atmodel Salbertea efairby 
good one. 

Densities for the input model were computed from 


(Birch, 1964) 


0 = 0.252+0.3788Vp 


and the shear wave velocities were calculated from the 


corresponding P-wave values by assuming the Poisson's ratio 
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computing the amplitudes. The synthetic obtained from the 
input model M1 from horizontal layer calculations is shown 
ieeeoLaure-6.o. LeeeiG  clearesthat. the tinsteanrivals are 
successfully reproduced by the synthetic. The Pn arrivals 
are very weak and the maximum amplitude for PmP lies at 140 
km which according to our previous calculation should be 132 
km. Also it does not reproduce all the reflection branches 
that we identified on the observed seismogram. Considering 
the results from the T’?-X? interpretation we then introduced 
two low velocity zones at 11-14 and 24-27 kilometers of 
depths. Most of the features were in good agreement except 
that the PmP branch was much earlier than is required by the 
data. Adjusting the low velocity zone at 24-27 km allowed 
US tGutix =thusevery easily we hiss led *tosamcritical distance 
closer to the shot point (figure-6.9) and was an additional 
advantage of the new model, since this agrees better with 
the observations. Actually the position of the maximum 
amplitude of the PmP branch depends on the frequency content 
of the corresponding signals. According to Cerveny (1966) 
PhemcEL EbtCcalspointe ton ay OeHertz Signal retlected  sinom the 
Moho interface usually lies at about 10-25 km beyond the 
theoretical value. Therefore observed and theoretical 
Chiprca eed 1 StanceSm Orem l ole andes 152 9kms “for the PmP are 
clearly admissible. The new model is shown in figure-6.10 
and will be called the model M4 hereafter. Table-6.4 
contains the adjusted parameters for the model M4. A ray 


diagram »corresponding to model M4 is shown in figure-G.11. 
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TABLE:6.4 Nominal parameters as calculated 
for the final model (M4). 
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Figure 6.8 Synthetic seismogram obtained from model M1 of 


horizontal layer interpretation (a dot on the PmP 


the position of the maximum amplitude). 
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Figure 6.9 Synthetic seismogram from the final model M4 (a 


dot on the PmP indicates the position of the maximum 


amplitude). 
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Figure 6.10 Final model (M4) showing the crustal “stucture 


under the study area. 
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Figure 6.11 Ray diagram corresponding 


tinal) model 
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Dashed curves represent the wave fronts. 
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We consider model M4 as our final presentation of the 
Crustal structure under the study area. Of course some of 
the features remain quite speculative. This is a limitation 
we could not avoid due to the lack of a reversed profile and 


the difficulties associated with obtaining amplitude 


information. 
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7. DISCUSSION 

In the last chapter interpretation of the observed data were 
Canpredsscut  oneethe basis of Vfinsteandireflectedtarrivais? 
Critical points and approximate relative amplitudes selene ws 
also worthwhile to discuss the derived crustal model in 
terms of the tectonic history of the region, the lithology 
at various depths and if possible the physical parameters 
such aS pressure, temperature, conductivity, porosity etc as 
eecunctionmsor depth. 

Variation of the seismic velocities with depth in the 
crust can be related to the changes in gross chemical 
composition and/or progressive metamorphism. Many 
geo-scientists consider that generally the earth's crust 
consists of an upper crystalline basement composed of 
metamorphic rocks such as gneisses and schists, an 
intermediate zone of migmatites and the lower crust of 
granulite facies rocks. The observed rapid velocity 
increase from 5.97 to 6.2 km/sec in the basement rocks may 
correspond to a gradient §zone with increasing closure of 
microfractures due to overburden pressure. Beneath this 
gradienteezone Stheavelocity becomes nearlyfoconstantlatmond? 
km/sec. This fact iS consistent with the behavior of either 
doyemorsusaturatedimrocks Wsubjectedstomincreasangrconfinang 
pressure, pore! pressure) being dow.) «The next’ layer. of » 6.7 
km/sec velocity is a Similar’ zone with a difference in 
chemical composition and/or gradational metamorphism. 


Crustal velocities of the order of 6.0 km/sec are 
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mmterpretedias an indicationoftqranitic rocks®¢while © ithose 
with@emore §than= 7.0°skm/sec “are vasstmed®to be ‘of “gabbroic 
composition. A velocity of 6.7 km/sec in the middle and 
VOVermmcvUStm@mthereforep ishould@icornespondm towearlayer of 
iNteEMedtale composi@uvens. The lowe evelocity min wathnemealower 
part of the basement was suggested by several factors: 
Gi@rapidwattenuation of thesPgq branchBafter. 70 ckm: 

(2) observed reflection branch from a possible interface at 
11 km depth. 

(3) observed high amplitude of the refracted/reflected 
Signals from the top of the 6.49 km/sec layer. 

(4) observed critical point much smaller than the calculated 
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ene homgche @ceblection branch iromeethe «Stoperot Pathe 
km/sec layer (theoretical and observed values are 109 & 90 
kms respectively). 

A low velocity zone at the lower basement may be 
explained as due to an intrusive sheet and possibly relates 
to a high conductive region in the crust as proposed by 
Gamerela, andeGough nel! SP. )er Fromefiqure-2'48 rt appears ¥ that 
the NACP anomally is located approximately between 
longitudes 104.5°W and 105.5°W at latitude 49°N. This means 
that normally we should expect to see any seismic expression 
Gambheeanomalblyiwethines0=10Sekm from the shot point. From 
the observed seismogram we find that a rapid attenuation of 
the Pg amplitudes occur between 65 and 100 kilometers’ from 
the shot. This suggests that there may be an inhomogeneous 


structure present in the crust ‘under the basement ©© From the 
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interpretation of these attenuated and delayed phases we 
proposed a low velocity structure under the 6.2 km/sec layer 
Gsechion Gal Amiliiwegwantetesrelate this velocity low -to 
the NACP anomally the following discrepancies arise: 
(1) In north Saskatchewan, the narrow conductive anomally 
was shown to be associated with steeply dipping graphite 
Sheets, which means a near vertical structure of the 
conductive body. The low velocity zone obtained from this 
Study was only a thin horizontally layered structure. 
(2) the NACP body was proposed to be at the surface of the 
Precambrian. The present study indicated a low velocity 
zone just below the basement layer. 
(3) suggested velocity for the low velocity zone was of the 
order of 6.1 km/sec from this study whereas the conductive 
materials of Camfield and Gough were thought to be graphitic 
sheets having a much lower velocity. 
However the following points are also important: 
(1) the low velocity zone was indicated by the anomally in 
Bhemdatanjust at the regqicn) (wheres =the | NACP “anomally was 
expected. 
e TetHemmeceDunmEOremcnew concictive Bbodygeuat Bagmpanuicutar 
location does not necessarily imply that it will have the 
Same depth at other locations. 

In fact the NACP body was. shown to have different 
depths between 3 and 75 km at different locations of North 
America. Therefore amesiugqhtiy eqreateraugdepuhh got (othe 


conductive body at the study area is clearly admissible and 
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we can conclude that the observed low velocity below the 6.2 
km/sec layer might correspond to the same conductive body as 
proposed by Camfield and Gough. 

The velocity inversion at 25-29 km depth is more well 
developed and is suggested by the computed interval 
VELOC RCVeS@EEECI@Et ck” Ranalysis ©) Thesvelocityedowsatethis 
depth could also be related to high conductivity as 
mentioned by Sternberg and Clay (1977), Barry et al (1977) 
and Chaipayungpun and Landisman (1977). The existence of a 
low velocity layer in the middle crust can be related to 
concentrated waters in the first formed granitic melts. it 
iS quite possible that this ‘wet zone' may have a velocity 
lower than the material above and even lower than _ the 
dehydrated material below (Berry, 1972). The release of 
free water into the surrounding pore spaces would provide a 
consistent explanation for the lowered seismic velocities, 
the increase of attenuation; the” reduction “of cohesive 
Strength and also the higher electrical conductivity in the 
CEUSE. 

FOrmESiMpliGity we@tooksatheserust -Mohomtransitioneto be 
Seuursteorder drsconminuiey, but@we didstnet. Yruve) Pour. the 
POsciDulTEye Comeame broader = stEalisiPiOnezonerat sthis#reqion® 
The observed characteristics such as long duration of the 
Signals, high amplitude reflect tons @arrivals® and =the 
cpl ieuingeo aabheeePmPesrerleécevons emays=Suggest™sa complex 
Paverwngqpecbesternate loweand highyvelocities (michs) 969; 


Clowes & Kanasewich, 1970; Giese 1976). Moreover numerical 
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experiments indicated that the only effective way to 
increase the apparent time duration of the reflected energy 
is to increase the thickness of the transition zone. 

Considerable effort was spent in attempting to find 
other models that would explain all the main features of the 
observed seismogram. All such efforts were unsuccessful. 
Therefore we Suppose that the model presented in this’ study 
1S a reasonable representation of the crust as illuminated 
by a source and viewed along a line in one direction from 
the source. The result could be more authentic if we had at 
our disposal, data from a reversed profile along the same 
line of observation. Refraction studies from adjacent 
profiles indicated that apaee from the local deviations the 
dip of the crustal layers hardly exceeds 3°, therefore true 
values of depths should not be significantly different from 
the calculated ones. In any event, the general structure of 
the layered crust should remain more or less the same. 

We also compared the derived final model with the 
Cestltseon OLhenerefraction™ prourles@icloseWatom thetistudy 
area. Comparing to the Swift Current and Weyburn models in 
figuses2 wimwe aSee athavrertiossupper Hcrustaleeveloci biess tare 
Similar. A refractor having a velocity of 7.1 km/sec and 
depuicmemecou sand: 4 kmssrespectively) “as “identified vat 
Suffield and Weyburn is not reflected in our data, instead 
we found a shallower layer with a velocity of 6.7 km/sec. 
Tt uomepossi bles thatmerheses twousvelocities Bhavessimilan 


Gompositionsaw sthedi fierence in densities. elt waseeproposed 
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that (Green et al, 1979) a crustal fault at 103°W is 
responsible for the differences in Moho depths on either 
Side of it. Therefore even without the knowledge of the 
correct position of the fault we could possibly explain the 
variation of the Moho depths at Weyburn and Limerick as due 
to this fault. The difference between Moho depths at Swift 
current and Limerick (about 2.0 km) is small and can be 
associated with a small dip. It again proves that the 
thickness of the crust does not vary too much at these two 
localities. 

Figures 7.1 & 7.2 show the interpreted results from the 
hol Seenonrthosouth and S197? Meast-westererracts on prolriies. 
The depths and velocities of the refractors obtained from 
our study 1S supported by these results. The results 
Obtainedetrom=mthe 1979 morth-southsprotiles are similar to 
the ones obtained from this study only for the upper and 
middle crust. An increase of moho depth by 3 km towards the 
south can be related to the Williston Basin tectonics. 
Comparison of the refraction results aS a whole indicated 
thaws 
(1) our profile lies entirely within the Churchill province. 
(2) the simple layered model of the crust is well suited for 
this area. 

(3) the structure iS more or less homogeneous in the 
vertical direction. 

The position of the 103°W fault could not be identified 


from our data since the profile does not really cover’ the 
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Figure 7.2 Crustal model from CO-CRUST 1977 east-west 


refraction profile (after Moon and de Landro, 1981). 
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fault area (figure-1.1). However combining the results of 
this study and the ones’ from the 1977 east-west profile 
\Mocumandscealandro, —19G ))minducatesethat “the, position sof 
the fault should be somewhat to the east of longitude 103°W. 
Determination of the exact position of the fault may be 
possible by a near vertical reflection experiment in an area 
between longitudes 103°W and 100°W near the international 
border. 

Due to inadequate control over the seismometer gains we 
had difficulties in interpreting the dynamic characteristics 
of the waves. Therefore we recommend that more care should 
be taken in future to retain the amplitude information. To 
achieve a reasonable control over the sedimentary velocity, 
more detectors should be placed within 15 kilometers of the 
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